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FIXED-BASE VISUAL SIMULATION OF OBSTACLE AVOIDANCE 
DURING TERMINAL DESCENT OF ADVANCED APOLLO SPACECRAFT 
WITH AN ALL-FLEXIBLE PARAWING 

By G. Kimball Miller, Jr. 

Langley Research Center 

SUMMARY 

A fixed -base visual simulation study using wind-tunnel aerodynamics has been con- 
ducted to determine the ability of an onboard pilot to control a single -keel all-flexible 
parawing and an advanced Apollo spacecraft combination so as to avoid obstacles in the 
desired landing area. In the investigation, a closed-circuit television system in conjunc- 
tion with an earth -terrain model was employed for image generation, and six rigid-body 
degrees of freedom of the vehicle were simulated. The pilot controlled the vehicle 
through a single-axis hand controller which activated constant-rate reels to induce dif- 
ferential changes in the wing -tip lines for lateral control. Longitudinal and lateral trim 
controls could be commanded through a thumb switch on the hand controller. The simu- 
lated viewing system provided the pilot with a field of view of 34.6° vertical by 48.3° hor- 
izontal. The pilot's task was to assume control of the vehicle under the influence of 
winds with maximum velocities up to 26 ft-sec“l (7.9 m-see - *) at initial altitudes up to 
2000 ft (609.6 m) and land as close as possible to the center of a specified landing site 
with the vehicle headed into the wind to reduce ground speed from 29 ft-sec - l 
(8.8 m-sec"l) experienced in still air. 

The results of the investigation showed that the pilot avoided the obstacles 95.5 per- 
cent of the flights and landed with a heading error of 10° or less 59 percent of the flights 
with a miss distance generally less than 500 ft (152.4 m). Increasing the vehicle turn- 
rate capability by a factor of about 3 had little effect on terminal accuracy and resulted 
in a tendency of the pilot to overcontrol when the higher turn rates were used. 

INTRODUCTION 

Spacecraft recovery systems which employ all -flexible parawings to provide accu- 
rate terminal control and low horizontal and vertical velocities at touchdown are currently 
of interest. The ability of an onboard pilot to control an all -flexible parawing and 
advanced Apollo spacecraft combination at initial altitudes up to 18 000 ft (5486.4 m) so 



as to attain specified landing sites was demonstrated in a simulation study reported in 
reference 1. The reference indicated that a maximum turn-rate capability of approxi- 
mately 9° sec~l was satisfactory and that the pilots could consistently attain the desired 
landing site under various wind conditions. Although these low turn rates are generally 
satisfactory in descending flight, they may not be satisfactory for low -altitude obstacle 
avoidance. 

The present fixed-base simulator study was performed to determine the ability of 
an onboard pilot to control an all-flexible parawing and advanced Apollo spacecraft com- 
bination so as to avoid obstacles with near-zero vertical dimension, such as rivers, 
gullies, and fences, which are representative of those present in contemplated landing 
areas for this type of descent system. The flights were initiated at an altitude of approx- 
imately 2000 ft (609.6 m) with wind velocities from 0 to 26 ft-sec - l (7.9 m-sec-1) and 
various wind directions. The investigation employed the 1962 U.S. Standard Atmosphere. 

Pilots seated in the Apollo spacecraft cannot see the ground when the vehicle is 
suspended with the heat shield down. Consequently, an onboard viewing system was 
included which had a field of view of 34.6° vertical by 48.3° horizontal and could be con- 
trolled in pitch by the pilot. The pilot had a view from the local vertical forward 34.6° 
when the viewing system was directed in its maximum downward position. In its maxi- 
mum upward position, the pilot had a view of ±17.3° about the earth's horizon. 

Equations of motion permitting six rigid-body degrees of freedom were used to 
represent the vehicle and were solved in real time by a digital computer. The pilot con- 
trolled the vehicle through a single -axis hand controller. The inputs from the hand con- 
troller actuated constant-rate reels which changed the right and left parawing-tip lines 
for lateral control. 


EQUATIONS OF MOTION 

Measurements for this investigation were made in U.S. Customary Units but are 
also given in the International System of Units (SI). (See ref. 2.) The symbols used in 
the present study are presented in appendix A. The equations of motion which permitted 
six rigid-body degrees of freedom of the vehicle and the transformation matrices for the 
assumed axis systems are presented in appendix B. The rigid-body assumptions were 
made on the basis of observations of wind-tunnel models and flight-test models which 
indicated that the parawing and the payload are tightly coupled in fully deployed gliding 
flight. The form of the equations of motion was dictated by the requirements of the sim- 
ulation equipment used in the study. The force equations were written with respect to 
cylindrical coordinates and the moment equations with respect to body axes. 
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To permit the simulation of several different configurations without extensive 
reprograming for each configuration, the computer program used in the present inves- 
tigation included the equations (appendix C) required to transfer the basic aerodynamic 
derivatives of the parawing from the confluence of the suspension lines to the center of 
gravity of the parawing -payload combination. The equations added to the flexibility of 
the program in that simulated configurations could be completely changed merely by 
changing a few constants. 


VEHICLE DESCRIPTION 

The simulated configuration consisted of a single -keel 45° all -flexible parawing 
with a one-eighth keel -length nose cut and an advanced Apollo spacecraft suspended as 
indicated in figure 1. The simulated all -flexible parawing was geometrically similar to 
that described in reference 3, including the number and placement of the suspension lines. 
The wing loading of the vehicle (parawing -payload combination) was 1.28 lb-ff 2 
(61.28 N-m“ 2 ). The wing material was assumed to have a weight of 2.2 oz-yd~2 
(0.075 kg-m-2) and the suspension lines had a weight of 128 lb (58.06 kg). Pertinent 
characteristics of the wing when flat were 


131.8 ft (40.2 m) 
b = 186.2 ft (56.7 m) 

S = 12 000 ft 2 (1114.8 m 2 ) 

The moments of inertia of the 15 000-lb (6804 kg) payload were 
*XP = slug-ft 2 (9897.3 kg-m^) 

Iyp = 6456 slug-ft^ (8753.0 kg-m^) 

*ZP = ^858 slug-ft 2 (7942.3 kg-m2) 

The cross products of inertia were assumed to be negligible and were omitted. 

The mass and inertia of a parawing is greater in air than in a vacuum because of 
the air trapped within the canopy envelope and because some of the surrounding air is 
also set in motion. The inertia of a 24-ft (7.3 m) keel parawing has been determined 
experimentally (ref. 4) and includes apparent mass and inertia effects. The inertias of 
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the wing and suspension lines, including apparent inertia effects, for the 131.8-ft 
(40.2 m) keel parawing used in the present simulation were estimated from the small 
scale data of reference 4 and are presented in figure 2. The apparent-mass effects were 
generally assumed to be zero in the present study. The parawing inertias (fig. 2) and the 
payload inertias were combined to give the total inertias of the vehicle by using the equa- 
tions presented in appendix C. 

During descent under the inflated parawing, the heat shield was down and the astro- 
naut's view through the advanced Apollo windows was upward in the direction of the wing. 
It was, therefore, necessary to include an auxiliary viewing system that would enable the 
pilot to view the earth. The viewing system used in the simulation had a field of view 
of 34.6° vertical by 48.3° horizontal. The direction of the optical axis of the viewing 
system could be controlled by the pilot to any angular position between the Apollo X-axis 
and a point approximately 75° below the X-axis. 

Because longitudinal control was quite small and thus of little practical use (ref. 1), 
only a single-axis hand controller was employed which provided lateral control by actu- 
ating simulated constant-rate ^0.02^ sec _ l) reels to deflect differentially the right and 
left wing -tip lines. A four-position thumb switch on the hand controller was used to trim 
the vehicle both longitudinally (by deflecting the rear keel line) and laterally (by differ- 
ential wing -tip deflection). 


SIMULATION EQUIPMENT 

The equipment used in the simulation is depicted in figure 3. A fixed-base cockpit 
(fig. 4) was fitted with a television monitor that provided a field of view of 34.6° vertical 
by 48.3° horizontal to the pilot. The instruments used during the flights were the three- 
axis gyro horizon, the altimeter, and the control line position indicators. The remaining 
instruments were used only for checkout purposes. A single -axis hand controller with a 
dead band that was 25 percent of maximum deflection in roll was located to the right of 
the pilot's seat. Depressing the trigger and deflecting the control stick activated 
constant-rate reels to drive the right and left wing-tip lines differentially. When the 
controller was returned to the neutral position with the trigger depressed, the control 
lines remained at their given position. The control lines, could be returned to the trim 
position by releasing the trigger with the control stick in any position. Lateral and longi- 
tudinal trim inputs activated the wing-tip-line reels and rear -keel -line reel, respectively, 
at O.OU^ sec - !. Control inputs were transmitted to a digital computer which solved the 
equations of motion in real time and generated drive signals for the simulation equipment. 
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The visual presentation was obtained by means of an 875-scan-line image orthicon 
television camera in conjunction with an optical pickup similar to that described in 
reference 5. The optical pickup was driven by the output of the moment equations to 
provide the three rotational degrees of freedom of the vehicle. The three translational 
degrees of freedom were obtained by mounting the optical pickup and camera combination 
on a transport system that moved relative to a terrain model in response to the output 
of the force equations. The earth -terrain model was a 22.0- by 34.9-ft (6.7 by 10.6 m) 
translucent back -lighted screen, scaled at 1200:1, and was an adaptation of part of a 
Langley Research Center simulator originally designed to study lunar operations. A 
transparent photograph of an artist's concept of the Wallops Island, Virginia, runway 
complex, representing an area approximately 12 000 ft (3657.6 m) square, was mounted 
on the front of the model (fig. 5). In addition, 20 appropriately scaled buildings were 
placed on the transparency. It should be noted that model protection considerations 
imposed a minimum altitude limit of approximately 30 ft (9.1 m), at which altitude an 
electrical stop was activated. 

AERODYNAMIC PARAMETERS 

The aerodynamic and control parameters employed in the present investigation are 
the same as those used in the study in reference 1 and do not include the payload aero- 
dynamics. All parameters are based on the span and area of the parawing when flat and 
are referenced to the system of reference axes (fig. 1) located at the theoretical conflu- 
ence of the suspension lines. 

The aerodynamic data were programed on a digital computer in the straight-line 
segments shown in figure 6. The computer program transferred the data to the center 
of gravity of the parawing -payload combination by means of the equations presented in 
appendix C. The data in figures 6(a) and 6(c) were obtained from the static wind-tunnel 
measurements of an 18 -ft (5.5 m) keel model presented in reference 3. 

The sideslip and dynamic derivatives in figure 6(b) were obtained from forced 
oscillation tests of a 5 -ft (1.5 m) keel model in the Langley full-scale tunnel. The 
small model was used because of its adaptability to existing test equipment. It should 
be noted that the oscillatory test technique used to obtain these data yields combination 
derivatives in which the acceleration terms appear. Thus, the derivatives in figure 6(b) 
are combination derivatives /of the form C, + C, . sin on and were obtained from mea- 

\ *P H ) 

surements of the force and moment components in phase and out of phase with the model 
motion. Since measurements of the acceleration derivatives do not exist at present, sep- 
aration of the combination derivatives was impossible. Reference 6, however, indicates 
that the aerodynamic effect of the acceleration derivatives for rigid leading-edge conical 
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parawings appears to be negligible in the angle -of -attack range below the stall. There- 
fore, the oscillation test data were assumed to be simply those noted in figure 6(b), that 

is, the combination C, + C, sin a was assumed to be C 7 , and so forth. 

h l p L V 

TASK DESCRIPTION 

The flights originated at altitudes up to 2000 ft (609.6 m) with the vehicle upwind 
of a designated landing site and generally headed into the wind. The pilot's task was to 
control the vehicle so as to attain the center of the designated landing site with the vehicle 
heading into the wind at flight termination (30-ft (9.1 m) altitude) to reduce ground speed. 
The primary landing site used in the present study was the 500-ft by 700-ft by 750-ft 
grassy triangle (152.4 m by 213.4 m by 228.6 m) formed by three intersecting runways 
(fig. 7). The minimum turn radius of the advanced Apollo and parawing combination was 
approximately 225 ft (68.6 m) (ref. 1). A circle of this minimum turn radius can not be 
accomplished within the boundary of the primary landing site (fig. 7). Because the vehicle 
was not capable of turning 360° inside the triangle, the acceptable landing area included 
the three grassy plots adjacent to the triangle. The runways separating these four grassy 
plots constituted obstacles to be avoided. These two-dimensional obstacles were believed 
to be representative because the landing areas contemplated for the advanced Apollo and 
parawing descent systems contain obstacles with a near-zero vertical dimension; that is, 
the landing area will probably contain gullies, streams, or fences, but will not contain 
forrests or high buildings near the primary landing site. 

Four wind profiles (fig. 8), scaled in decreasing magnitude from A to D, were used 
in the investigation. The wind magnitudes ranged from very nearly zero for wind D at 
low altitude to approximately the still air forward speed of the vehicle (29 ft-sec-1 
(8.8 m-sec-1)) for wind A at flight initiation. For a given wind profile, the flight initia- 
tion point was always chosen so that the primary landing site could be attained with proper 
maneuvering. In all flights, the terminal magnitude and direction of the wind were known 
by the pilot. 


RESULTS AND DISCUSSION 

The results presented herein were obtained by a pilot (research engineer) who was 
considered to be well trained inasmuch as he had participated in the simulation study 
reported in reference 1 and in several radio-controlled flight tests of small-scale 
parawings. 


6 



Preliminary Results 


Although 131.8-ft (40.2 m) keel parawings have not been flight tested at the present 
time, a 24-ft (7.3 m) keel parawing has undergone preliminary flight tests. Therefore, 
to provide a ; comparison between simulator results and flight -test results, a parawing 
with a 24-ft -long (7.3 m) keel and a wing loading of 1.15 lb-ft“2 (55.06 N-m“2) was sim- 
ulated by using the aerodynamic parameters presented in figure 6. At the present time 
flight -test data at a wing loading of 1.15 lb~ft“2 (55.06 N-m-2) are not available. Con- 
sequently, the flight -test data presented herein for a wing loading of 1.15 lb-ft -2 
(55.06 N-m-2) were interpolated from flight -test data obtained at the Langley Research 
Center at wing loadings of 0.394, 0.700, 1.000, and 1.300 lb-ft"2 (18.86, 33.51, 47.87, and 
62.24 N-m-2). There is good agreement (fig. 9(a)) between turn rates (obtained by timing 
360° turns) experienced during flight tests and those generated on the simulator for lat- 
eral control inputs up to 0.033f k , the maximum input employed with the small-scale para- 
wing. In addition, the increase in rate of descent with an increase in turn rate obtained 
by simulation (fig. 9(b)) agrees reasonably well with the flight -test data. It thus appears 
that the aerodynamic parameters used in the simulation are representative of those for 
all -flexible parawings. 


Apparent Mass and Inertia Effects 

In addition to the 131.8-ft (40.2 m) keel parawing with moments of inertia given in 
figure 2 and zero apparent mass, two other configurations were simulated to examine the 
effects of apparent mass and inertia variations on pilot performance. The second con- 
figuration arbitrarily included an apparent mass of approximately 250 slugs (3648.5 kg) 
with no inertia changes, whereas the third configuration employed zero apparent mass 
and moments of inertia that were reduced from those of figure 2 by approximately 
10 percent for I^^y, 20 percent for Iyw> anc * ^ P ercen t for *zw‘ Time histories of 
a specific maneuver with each of the three configurations are shown in figure 10. The 
pilot was unable to detect the differences in the flight characteristics of the three config- 
urations in figure 10. However, the vehicle turn rate for a given lateral control input 
was somewhat less for these three configurations than it was for the vehicle with zero 
apparent mass and inertia reported in reference 1. For example, a control input 
ALq/Zjj of 0.04 resulted in an approximately 4° sec -1 turn rate for the three configu- 
rations of figure 10, whereas the vehicle in reference 1 was capable of a 5.5° sec~l 
turn rate when the apparent mass and inertia were omitted. 

Piloting Techniques 

Ground tracks of representative flights made under the influence of the four wind 
profiles used during the study are presented in figure 11. The flight techniques used 
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by the pilot were generally the same as those given in reference 1. In the higher velocity 
winds A, B, and C, the pilot generally controlled the vehicle to some point downwind of 
the landing site, turned into the wind to land, and performed S-turns to lose altitude if it 
appeared the vehicle was going to overshoot the center of the landing site. The use of 
S-turns to lose altitude near the landing site was ineffective in low winds because the 
ground speed was relatively high and resulted in the vehicle moving farther away from 
the landing site than the pilot desired. Consequently, for wind D the pilot generally 
remained outside the boundaries of the primary landing site; he kept the landing site in 
view until altitude was reduced so that he could turn and run straight toward the. center 
of the target for touchdown. Occasionally, during the final approach, the pilot could not 
attain the primary landing site and it was necessary to land in one of the three alternate 
areas. Time histories of typical flights made in a low wind (wind D) and a high wind 
(wind A) are shown in figures 12 and 13, respectively. 

Terminal Accuracy 

The primary concern of the present investigation was to determine the ability of a 
pilot to avoid obstacles that might be located in the desired landing area. Approximately 
200 flights were made under the influence of various wind conditions (fig. 8). The pilot 
was able to avoid the obstacles (runways) 95.5 percent of the flights (table I) and had a 
terminal heading error of less than 10° (which was believed to be a reasonable allowable 
error) approximately 59.0 percent of the flights. The pilot performed emergency landings 
in one of the three alternate sites 11.4 percent of the flights. 

TABLE I.- TERMINAL ACCURACY 


Wind 

profile 

Number of 

Miss distance, ft 

(m) 

Heading error, AA, deg 

Percent of runs with 
heading error S 10° 

Percent of runs 

simulator 

runs 

Arithmetic 

mean 

Standard 

deviation 

Arithmetic 

mean 

Standard 

deviation 

Avoided 

obstacles 

Emergency 

alternate 

A. 

56 

191.7 

(58.4) 

164.4 

(50.1) 

-3.3 

19.7 

51.8 

94.7 

8.9. 

B 

58 

194.6 

(59.3) 

145.5 

(44.3) 

-1.2 

18.0 

60.3 

94.8 

13.8 

C 

31 

169.7 

(51.7) 

155.2 

(47.3) 

-4.8 

22.4 

71.0 

100.0 

6.4 

D 

56 

338.4 

(103.1) 

205.9 

(62.7) 

Not 

applicable 

Not 

applicable 

Not 

applicable 

! 96.5 

16.1 

All flights 

201 

229.9 

(70.1) 

182,2 

(55.5) 

-2.9 

20.3 

59.0 

95.5 

11.4 


Because little reduction in ground speed could be achieved by heading into the wind 
for low wind D, the pilot ignored his heading angle at flight termination and concentrated 
on landing as close as possible to the center of the landing site. Even though the pilot 
was able to concentrate more heavily on attaining small miss distances, the terminal 
accuracy experienced with wind D was considerably poorer than that with higher winds 
because the relatively high ground speed resulted in a tendency to overshoot the desired 
landing site. 
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In this investigation, the miss distances were considerably smaller than those of 
the simulation in reference 1 for which the mean miss distance was 618.6 ft (188.5 m) 
with a standard deviation of 446.4 ft (136.1 m). The landing area used in reference 1 had 
no recognizable landmarks near the primary landing site and the pilots frequently became 
disoriented at low altitudes. In the present study the pilot was able to compensate for the 
small field of view and obtain smaller miss distances by using the known pattern of the 
runway area as orientation cues relative to the primary landing site. During low -altitude 
maneuvers, the primary landing site was outside the small -angle field of view for extended 
periods of time. However, with a well-defined pattern of readily recognizable landmarks 
surrounding the primary landing site, the pilot was always able to remain properly 
oriented. 


Supplementary Investigation 

It was believed that a vehicle with a high -turn -rate capability (that is, which could 
turn inside the boundaries of the primary landing site) would make the obstacle avoidance 
task considerably easier. A turn-rate capability of approximately 25° sec"* was 
obtained by simulating the 24 -ft (7.3 m) keel parawing and using the aerodynamics param- 
eters presented in figure 6. 


TABLE H.- TERMINAL ACCURACY (HIGH-TURN-RATE CAPABILITY) 


Wind 

profile 

Number of 

Miss distance, ft 

(m) 

Heading error, AX, deg 

Percent of runs with 
heading error ^ 10° 

Percent of runs 

simulator 

runs 

Arithmetic 

mean 

Standard 

deviation 

Arithmetic 

mean 

Standard 

deviation 

Avoided 

obstacles 

Emergency 

alternate 

A 

40 

160.9 

(49.04) 

126.5 

(38.6) 

-0.6 

31.9 

47.5 

100.0 

7.5 

B 

j 

40 

142.8 

(43.5) 

131.0 

(39.9) 

-8.7 

25.3 

55.0 

92.5 

7.5 

C 

41 j 

131.7 

(40.1) 

112.1 

(34.2) 

.9 

18.8 

65.8 

95.1 

4.9 

D 

40 

233.6 

(71.2) 

146.6 

(44.7) 

Not 

applicable 

Not 

applicable 

Not 

applicable 

100.0 

1 2.5 

All flights 

161 

167.0 

(50.9) 

133.9 

(40.8) 

-2.7 

26.0 

57.9 

96.9 

5.0 


The terminal accuracy that the pilot was able to attain with the high -turn -rate capa- 
bility is shown in table n for the four wind profiles used in the present study. The pilot's 
ability to avoid obstacles with the vehicle having a high -turn -rate capability was only 
slightly improved (96.9 percent successful) over that with the full-scale vehicle having a 
low-turn-rate capability (95.5 percent successful). It was not necessary, however, for 
the pilot to perform emergency landings in the alternate sites as frequently with the high- 
turn-rate vehicle, which resulted in smaller miss distances particularly for flights made 
under the influence of low wind D. The terminal heading error increased somewhat with 
high turn rates, primarily because of a tendency to overcontrol. To minimize the over- 
control problem that existed when high turn rates were employed, the pilot generally 
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planned the flights with turn rates of approximately 8° sec 'I and used the higher turn 
rates only in the event of an emergency. 

CONCLUDING REMARKS 

A fixed-base visual simulation study has been conducted to determine the ability of 
an onboard pilot to control a 131.8-ft (40.2 m) keel all -flexible parawing and advanced 
Apollo spacecraft combination so as to avoid obstacles located in the desired landing 
area. The investigation included six rigid-body degrees of freedom of the vehicle. The 
pilot's task was to control the vehicle, under the influence of winds with maximum veloc- 
ities up to 26 ft-sec-1 (7.9 m-sec-1), at initial altitudes up to 2000 ft (609.6 m) and land 
as close as possible to the center of a specified landing site with the vehicle headed into 
the wind to reduce ground speed. 

Within the limits of the simulation, the results of the study showed that the pilot 
could avoid obstacles with near-zero vertical dimension, such as rivers, gullies, and 
fences, 95.5 percent of the flights and land with a heading error of 10° or less approxi- 
mately 59 percent of the flights. A landing area that included a readily recognizable pat- 
tern of landmarks in the vicinity of the primary landing site provided orientation cues 
which were useful during low -altitude maneuvers when the primary landing site was out- 
side the pilot's narrow field of view. These orientation cues permitted the pilot to reduce 
the miss distances to less than one-half those experienced in the investigation of NASA 
TN D-5049 for which orientation cues were not available. 

To compare the simulated and flight -test vehicle characteristics, a 24 -ft (7.3 m) 
keel vehicle was also simulated. The turn-rate capability of the small-scale vehicle was 
approximately three times that of the full-scale vehicle and the simulated turn rates 
agreed very well with flight-test results up to about 25° sec - l. However, very little 
improvement in the pilot's ability to avoid obstacles and attain small miss distances was 
experienced by using high turn rates. The pilot generally planned the flights with turn 
rates of about 8° sec _ l and used the higher turn rates only in the event of an emergency. 
The pilot was unable to detect changes in vehicle performance due to inertia variations of 
approximately 10, 20, and 75 percent in roll, pitch, and yaw, respectively. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 30, 1970. 
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APPENDIX A 


"ij 


SYMBOLS 


The assumed axis systems are shown in figure 14. 


direction cosines, where i = 1, 2, 3 and j = 1, 2, 


span of flexible parawing when flat, feet (meters) 




Rolling moment about Xp-axis 
QSb 


- 


C n — 


C X - 


c y ~ 


c z = 


Pitching moment about Ypj-axis 
QSb 

Yawing moment about Zj^-axis 
QSb “ 

Force along X^-axis 
QS 

Force along Yj^-axis 


QS 

Force along Z^-axis 


QS 


C m 0 va l ue °f C m at zero angle of attack 
C x o va * ue °f C x at zero angle of attack 
Cz o va -l ue °f at zero angle of attack 


dC 


m 


'm 


a da 


■, per radian 


9C 


X 


X a da 


, per radian 



APPENDIX A — Continued 
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2V 


R 


9C n 

C n = — , per radian 

P 9 J2L. 

2V R 


9C, 


' Y P 9 JjbJ 


•, per radian 


2V 


R 


9C 


m 


'HI, 


* « qb 


per radian 


2V R 


9C V 
c = - A 
X q 9 Jib-’ 

2V R 



9C Z 

9_qb_ : 

2Vr 


per radian 


per radian 
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APPENDIX A - Continued 


9C 


C, = r — , per radian 

I-* „ rb 


^ 9 


2 V 


R 


9C n 

C n > per radian 

r a rb 


2V 


R 


9C, 


' Y r " 8 


per radian 


2V 


R 


9°X 9C Z 8C m 


AL 

7 

B 

AL 

7 

B 

AL 



l k 




variation in C^, C^,, and C m , respectively, due to lateral 
inputs 


9C V 9C n 9C 7 

A L A. 


9 


AL 

7 

B 

AL 

7 

B 

AL 

l k 




^k 


variation in C v , C n , and C, , respectively, due to lateral 

/3 /3 L p 

inputs 


9C X 

ac z 

ac 

m 

8^.’ 



*k 

Ik 

^k 

8C y 

9C n 

8C I 

9 4L’ 

8^’ 

8^L 

^k 

*k 

^k 


trim control parameters (that is, changes in coefficients due to 
longitudinal trim inputs) 


lateral control parameters (that is, changes in coefficients due 
to lateral inputs) 


c, d, e, f 
F XB’ F YB’ f zb 


quaternions or Euler parameters (ref. 7) 

aerodynamic force along Xg-, Yg-, and Zg-axis, respectively, 
pounds (newtons) 


S e 


acceleration at surface of earth due to gravitational attraction, 
32.2 feet-second - ^ (9.81 meters-second - 2) 
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APPENDIX A - Continued 


altitude above surface of earth, feet (meters) 


X XB’ ^TB’ *ZB 


X XP’ X YP’ *ZP 


moments of inertia of parawing -payload combination including apparent 
inertias, slug-feet^ (kilogram-meters^) 

moments of inertia of payload, slug-feet^ (kilogram -meter s^ ) 


^XW’ *ZW moments of inertia of the parawing and suspension lines including 

apparent inertias, slug-feet^ (kilogram-meters^) 


AL 


sum of lateral control and trim inputs (that is, differential deflection of right 
and left wing -tip lines where right roll is obtained by shortening right and 
lengthening left tip lines (positive AL)), feet (meters) 


AL/ 


lateral control input, feet (meters) 


AL n 


lateral trim input, feet (meters) 


parawing keel length when flat, feet (meters) 


Al 


longitudinal trim input (that is, change in rear -keel -line length from that 
required for trim at a = 30°), feet (meters) 


M XB’ M YB’ M ZB aerodynamic moment about X B-> Y B-’ and Zg-axis, respectively, 

foot-pounds (meter -newtons) 

m a apparent mass, slugs (kilograms) 

mg mass of parawing -payload combination excluding apparent mass, slugs 

(kilograms) 

mp mass of payload, slugs (kilograms) 

mass of parawing and suspension lines excluding apparent mass, slugs 
(kilograms) 


p, q, r vehicle angular velocity about Xg-, Yg-, and Zg-axis, respectively, 
radians/ second or degrees/ second 
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APPENDIX A - Continued 

Q dynamic pressure, pounds -foot (newtons -meter -2) 

R, \f r , Y cylindrical coordinate system with origin at center of earth and vector R 
and angle SI' in XZ-plane (see fig. 14) 

R e assumed radius of earth, 20.9 x 10® feet (6.37 x 10® meters) 

S area of parawing when flat, feet^ (meters2) 

t time, seconds 

Vp relative wind velocity, feet- second - (meters-second - *) 

V XY velocity in plane parallel to surface of earth (ground speed), feet-second - l 

(meters - second - 1) 

X, Y, Z fixed-reference axis system with origin at center of earth which was 
assumed to be a nonrotating homogeneous sphere (see fig. 14) 

X B , Yg, Z B vehicle body axis system with origin at center of gravity of parawing- 
payload combination (see fig. 14) 

X c , Y^,, Zq moving-reference axis system with origin at surface of earth and 
with Zq -axis alined with local vertical and positive inward, 

X^-axis positive westward, and Y^-axis positive northward 
(see fig. 14) 

x o> y O’ z o optical axis system with origin at center of gravity of vehicle (differs 
from body axes by angle 6^ (see fig. 14) 

Xp, Yp, Zp axis system with origin at center of gravity of payload and parallel to 
vehicle body axis system 

X R’ Y R’ Z R orthogonal reference axis system with origin at theoretical confluence of 
parawing suspension lines and with Z R -axis alined with theoretical 
60-percent keel line (see fig. 1) 

X W’ Y W’ Z W ax * s system with origin at center of gravity of parawing and suspension 
lines and parallel to vehicle body axis system 

15 
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APPENDIX A - Continued 

y displacement along Y-axis, feet (meters) 

a angle between Xg-axis and projection of relative wind vector in XgZg -plane 

(referred to as angle of attack), radians or degrees 

/3 angle between relative wind vector and X B Z B -plane (referred to as sideslip 

angle), radians or degrees 

y flight -path angle (that is, angle between inertial velocity vector and 

XqY c -plane), radians or degrees 

7]p, ? 7 W , ? 7 b displacement along Z R -axis of center of gravity of payload, of parawing 
(including suspension lines), and of total composite body, respectively, 
feet (meters) (see fig. 1) 


9q angle of optical axis of viewing system in X R Zg -plane measured from 

Xg-axis, degrees 


X 


AX 


direction of flight referenced to north, defined as angle between 
and Y^-axis, degrees 


V X Y“ axis 


heading error (that is, difference between heading into wind and actual flight- 
path direction), degrees 


|p, £ w , |g displacement along X R -axis of center of gravity of payload, of parawing 
(including suspension lines), and of total composite body, respectively, 
feet (meters) (see fig. 1) 

p air density, slugs-foot"^ (kilograms-meter - ^) 

$ angle between Yg-axis and projection of relative wind vector in Y B Z B -plane 

(referred to as bank angle), radians or degrees 


ip, 6, <p Euler angles of rotation relating body axes and fixed-reference axis system, 
radians or degrees 

^X’ ^Y> wind component along X-, Y-, and Z-axLs, respectively, feet-second - * 

(meters-second -1 ) (see fig. 14) 
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APPENDIX A - Concluded 


Notation: 


absolute value 



square matrix 
column matrix 
transpose of matrix 



m=B,C,I,0 

n=B,C,I,0 


matrix which transforms a vector from axis system m to axis 
system n; B, C, I, and O represent body system, moving- 
reference system, fixed-reference system, and optics system, 
respectively 


A dot over a symbol indicates a time derivative. 

A prime denotes that the aerodynamic coefficient has been transferred to the center 
of gravity of the parawing -payload combination. 
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APPENDIX B 


EQUATIONS OF MOTION 

The simulation equipment used in the present study was a modification of the modi- 
fication of the lunar orbit and landing approach simulator at the Langley Research Center. 
The form of the equations of motion was a result of the requirements of this equipment. 

Force Equations 

The equations of motion for the three translational degrees of freedom in the cylin- 
drical coordinate system (fig. 14) are as follows: 


R - R'F^ = - 


fzc . /?e 

&e V R 


(Bl) 


RSF + 2R\& = 


XC 


m 


B 


(B2) 


Y = 


YC 


m 


B 


(B3) 


where F^q, FyC’ anc * ^ZC are aer °dy na mic forces along the X^-, Yq-, and 
Z^-axes and are given by 


with 


r ^ 

F xc 


f 

f xb 

\ f yc 

> - [ r c- B ] T < 

f yb \ 

F zc 

V. J 


f zb 
^ J 


F YTJ = QS|C' n + C’ a+ C’ + 

XB ^ I X,0 X a X q2V R 8 A i l k 

\ l k 


dC k 


AL 


AL 




(B4) 


(B5) 


18 



APPENDIX B - Continued 


F yb = Q s [ c y 


P 


+ 


Pb c . rb 
2V R Y r2VR 


8C y al 


9C y, 


|AL i k 


T~ ■* iTT7 P 


AL 


AL 

Ik 


(B6) 


and 


9C' 


F ZB = Q s fC’ o + C’ ce+C' -^-+— + 
ZB Z,0 Z a Zq 2V r 9 A l h 


9C Z 
9 AL 
*k 


AL 


(B7) 


The matrix relating the moving -reference axis system to the body axis system is orthog- 
onal and is 


where 



and 



The direction cosines are given 


Fc-b] " Fi-b] fc-l] 


(B8) 

cos ^ 0 

-sin ^ 


0 1 

0 

(B9) 

sin ^ 0 

cos ’3> 


A 11 A 12 

A 13 


A 21 A 22 

A 23 

(BIO) 

A 31 A 32 

A 33 


i terms of quaternions by 


A n = 2(c 2 + f 2 ) - 1 


(Bll) 

A 12 = ^( cc * + e ^) 


(B12) 

Ajg = 2(ce - df) 


(B13) 
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APPENDIX B — Continued 


A 21 = 2(cd - ef) 

(B14) 

A 22 = 2 (d 2 + f 2 ) - 1 

(B15) 

Ags = 2 (cf + de) 

(B16) 

A 31 = 2 (ce + df) 

(B17) 

A 32 = 2 (de - cf) 

(B18) 

A .33 = 2 (e 2 + f 2 ) - 1 

(B19) 

and 


c = ~ (fp - eq + dr) + Kgc 

Li 

(B20) 

• 1 

d = - (ep + fq - cr) + Kgd 

Li 

(B21) 

e = ~ (-dp + cq + fr) + Kge 

Li 

(B22) 

f = - — (cp + dq + er) + Kgf 
2 

(B23) 

where 


g = 1 - (c 2 + d 2 + e 2 + f 2 ) 

(B24) 

and K is a gain factor determined empirically on the computer. 

Moment Equations 


The equations of motion for the three rotational degrees of freedom in the body- 

axis system are as follows: 


P = + ( r YB ” X ZB) qr 

(B25) 



q = £^|^YB + ( r ZB ‘ I Xb) 1 'p] 

(B26) 
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APPENDIX B - Continued 


r “ iJof^ZB + ( ! XB " X Yb) P< ^1 
ZB — 


(B27) 


The aerodynamic moments are given by 


8C ? 3C j 

m yr = Q sb fc; p+ c: -sL+c; ^-+ — + _£_« 

XB I l B Z p2V R i r 2V R 9 ALi k + 8 |AL|^ 


7^11 (B28) 

l k 


sc* AJ ac: 


M vr = QSb/C^ 0 + C’ «+ C' -t+— ^ + 
YB | m,U m a m q 2V R d AL l k 


m 


AL 




AL 


(B29) 


and 


M 7R = QSbfCn /3+ Cl ^-+0^ 

ZB | n p H n p 2Vr n r 2y 


v. 9C' AT 

rb n AL 


8C* 


, AL 


A 


R 8 -k 8 
z k 


AL 


^k 


/3 


AL 


(B30) 


Auxiliary Equations 

The matrix which transforms a vector from the moving-reference axis system to 
the optical axis system is 


where 


[ r c-o] 


[ r B-°] = 


■sin 0 r 


1 

0 


ni 

*12 

ns 



*21 

r 22 

r 23 


(B31) 

I 31 

r 32 

CO 1 
CO 



cos Oq 

0 

sin e Q ~\ 



COS Or 


(B32) 
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APPENDIX B - Continued 


and are direction cosines relating the optical axis system to the moving-reference 
axis system. Thus, the angles i//, 9, and cp used to drive the optical pickup are 


sin 9 = -yjg 

+ J 1 ' ( y 13) : 


COS 


• t y 23 
sin cp = 


cos 4> = 


cos 9 

y 33 
cos 9 


Yio 

sin \p= 14 


cos 9 


- y ll 

COS !// = 


COS 9 


Dynamic pressure is given by 


Q = |py2 



(B33) 


(B34) 


(B35) 


(B36) 


where 


and 


V R ~ /( u + Scb) 2 + ( v + *\b) 2 + ( w + °ZBy 


(B37) 


n XB 


^X 


fi YB/ " (Ji-b) S / 


(B38) 


n ZB 
V. J 
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APPENDIX B - Continued 


r > 
U 


f. ^ 

s c 

V 

(“Mi 

& 

II 

Y c> 

W 
V > 


J 


with 

X c = Ri 

• • i 

Y C = Y / 

Z c = -R 


(B39) 


(B40) 


Angle of attack, sideslip angle, flight -path angle, bank angle, and altitude are 
expressed as follows: 


sin a = 


cos ci — 


sin /3 = 


cos /3 = 



w + 

°ZB 


vt u + 

°xb) 2 

+ (w + 

°zb) 2 


U + 

^XB 


\J(- 

"xb) 2 

+ ^w + 

Sb) 2 

v+ n YB 



V R 




\/( u + 

Scb) 2 

+ (w + 

"xb) 2 


•\ 




J 




J 


(B41) 


(B42) 
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APPENDIX B - Concluded 



(B43) 


(B44) 


(B45) 
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APPENDIX C 


TRANSFER EQUATIONS 

The computer program used in the present investigation included the following 
equations to transfer the basic aerodynamic derivatives, which were referenced to the 
theoretical confluence of the suspension lines (fig. 1), to the center of gravity of the 
parawing -payload combination: 

C' Y = C v (Cl) 

P */3 


i 

c v = C v - 2C V -2 (C2) 

Y P Y P Y /3 b 


, « B 
C v =C V +2C V -2 

*r *r y /3 b 

(C3) 

c = c - c S. 

*0 *0 /3 b 

(C4) 

, /^1r\ 

C, = C, - 2C, -2+2C v -2 

^p ^p J j3 b Y p b Y ^\b / 

(C5) 

C - C, + 2C, 13 C v 13 - 2C V a 13 

t-r ^ b Y r b b 2 

(C6) 

- ?R 

C n = C n + C v -2 
0 n /3 Y 0 b 

(C7) 

T T1t> Itd ^R^R 

C n P Cn p ■ 2Cn 0 b + C Yp b “ 2C Y /3 b 2 

(C8) 

c; r = C nr + 2C n ^ + C Yr ^ + 2C Y(3 (^) 

(C9) 

C X,0 + C X a a= C X,0 + C X a a 

(CIO) 
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APPENDIX C - Continued 


C’ = - 2C X 

A q A q A a b 


c 'z,o + c z a a = C Z,0 + C Z, a 


t ^ B 

C 7 = C z - 2C Z -f 

Z q Z q Z « b 


a 


ft + Cjjj 01 — C w ft + C m Q! + 


'm,0 


a 


m,0 


a 


( c x,o + c x a °j b “ ( c z,o + c z a t, B 


t ^?T5 

Cm = Cm - 2cm ~r — I- C-y — 2C 


€t 


a 


B 


' m q " m q 


y m 


a b T ^X q b 


is^B.c ^ +2C uj, 
x a b 2 z q b Z Q!\ b 


(Cll) 
(Cl 2) 
(C13) 
(C14) 
(C15) 


9C k 


8C-, 


AL 

3 

AL 

*k 


*k 


(C16) 


9C z 9C z 


AL 


AL 


(C17) 


9C m 9C m 


9C X 


8C 


Z ? B 


AL 


AL 


AL 


AL 


(CIS) 


9C X 9C X 

8 — 
l k Ik 


(C19) 


8C z _ 9C Z 

0^ 9 


(C20) 


8C m _ 8C m i 8C x % _ 9C Z % 
^T + b " 8 AL b 
I k z k 


(C21) 
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APPENDIX C - Continued 
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p, AL 

lk 


9 4L 
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dC i 

9c i 

3C y 
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l k 

AL 
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AL 
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8C' 

1 O 

9C V 
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AL 


AL 
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8C 


n. 
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AL 
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AL 
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Ik 

dC 7 
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AL 
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AL 

0 

AL 

9 

AL 

b 

*k 


lk 


^k 



(C22) 


(C23) 


(C24) 


(C25) 


(C26) 


(C27) 


where 


_/ m w V + f m P 

’ B \ m W + m p/ w \m w + nip 


m 


v B = 


w 


+ 


m w + m p / W 


f "P V 

im w + m p i P 


(C28) 

(C29) 


In addition, the moments of inertia of the parawing -payload combination about the 
Xg-, Yg-, and Zg-axes are given by 


T XB " X XP + m p(^B " ^p) 2 + r XW + ( m W + m a) (^W " ^b) 2 


(C30) 
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APPENDIX C - Concluded 
r YB = J YP + “P^B " ^p) 2 + (^B “ £p)^j 

+ *YW + ( m W + m a) [(% " ^b) 2 + (^B ‘ ^w) 2 ] 

X ZB = *ZP + m p(^B “ ^p) 2 + *ZW + ( m a + m p) (% " ^w) 2 


(C31) 

(C32) 
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3. Camera transport system 


Figure 3.- Simulation equipment. 
















Angle of attack, a, deg 


Angle of attack, a, deg 


Angle of attack, a, deg 


(a) Variation of longitudinal aerodynamic parameters and lift-drag ratio with a. 

Figure 6.- Aerodynamic and control parameters used in simulation. Data are referenced to axis system located at the theoretical confluence of the suspension lines as 

shown in figure 1. All parameters are based on area and span of parawing when flat. 
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Angle of attack, a, deg Angle of attack, a, deg Angle of attack, a, deg Angle of attack, a, deg 

(b) Variation of p, q, r , and p derivatives with a for undeflected controls. (These parameters are combination derivatives obtained from measurements made during 

forced oscillation tests of a 5-ft (1.52-m) keel model.) 


Figure 6.- Continued. 
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(c) Variation of control parameters with a. 
Figure 6.- Concluded. 











Turn rate, deg/sec 


(b) Rate of descent as a function of turn rate 


Fiaure 9.- Concluded. 
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(a) Wing with zero apparent mass and 
moments of inertia as shown in 
figure 2. 


(b) Wing with apparent mass of 
250 slugs (3648.5 kg) and 
moments of inertia as shown 
in figure 2. 


(c) Wind with zero apparent mass 
and reduced moments of inertia. 


Figure 10.- Apparent mass and inertia effects on vehicle flight characteristics. 


42 



Wind A 

Wind B 

Wind C 

Wind D 


Flight 0 
Initiation 
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Termination 



Figure 11.- Ground tracks of representative flights made under the influence of the four wind profiles used in the study. (Arrows denote 

wind direction at flight initiation.) 
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